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ABSTRACT 

Context. High-resolution X-ray observations in the imaging and spectral domain have recently opened a new window 
on active galactic nuclei (AGN) feedback onto the circumnuclear gas. Spectral diagnostics, as well as the remarkable 
morphological coincidence between [O ill] and X-rays, point to AGN photoionisation as the dominant ionisation mech- 
anism on scales as large as a few kpc. 

Aims. In this paper we extend these studies to the nearby Seyfert 2 galaxy NGC 1365, known to host a circumnuclear 
ring of intense star formation at ~1.3 kpc from the nucleus. The main scope of this investigation is to study the con- 
nection between nuclear activity and star formation in nearby AGN. 

Methods. We present a deep (~5.8 days) 0.3-2 keV high-resolution spectrum of NGC 1365, collected with the reflection 
grating spectrometer (RGS) on board XMM-Newton. 

Results. The spectrum is dominated by strong recombination lines of He- and H-like transitions from carbon to silicon, 
as well as by L transitions from Fexvu. The continuum is strong, especially in the 10 to 20A range. Formal fits require 
two optically thin, collisionally ionised plasma components, with temperatures ~300 and ~640 eV. However, they leave 
the bulk of the forbidden components of the He-a Ovu and Nvi triplets unaccounted for. These features can be ex- 
plained as being produced by photoionised gas. NGC 1365 is therefore the first obscured AGN, whose high-resolution 
X-ray spectrum requires both collisional ionisation and photoionisation. 

Conclusions. The relative weakness of photoionisation does not stem from the intrinsic weakness of its AGN, whose 
X-ray luminosity is ^ 10 42 erg s _1 . We suggest that it may instead come from the line-of-sight from the active nucleus 
to the NLR being blocked by optically thick matter in the broad line region, at the same time responsible for the large 
observed variation of the column density obscuring the X-ray active nucleus. Alternatively, NGC 1365 could host a 
remarkably luminous nuclear starburst when compared to the AGN accretion power. 

Key words. Galaxies: active - Galaxies:Seyferts - Galaxies: starbursts - X-rays:galaxies - X-rays:individual:NGC 1365 
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^ , 1. Ionising the NLRs: the X-ray view IDas et al. 20061 for a different view) raise the possibility 

- - that the interaction between radio jets and the NLR 

Most of the observational constraints on the physics and gas may affect the ionisa tion state of the NLRs as well. 

geometry of narrow line regions (NLRs) in active galactic Experimental support to t his hypothesis ha s so far been 

nuclei (AGN) come from optical spectroscopy and imag- presented in one case only IjAxon et al. 1998|> . 

ing. Comparison between spatially resolved spectroscopy 

of nearb y bright obscured AGN a n d photoionisation r n tnis cont ext, X-rays have recently opened a comple- 

models (|Ferland fc Osterbrock 19861 IQste rbrock 1989|) me ntary window on NLR physics. Whenever (extended) 

and the study of ionis ation cones ( |Pogge 1988[ NLRs in nearby Seyfert 2 galaxies are larger than the 

ITadhunter k Tsvetanov 1989|) indicate that the AGN angu lar resolution of the Chandra optics (~0.5"), X- 

high-energy emission is the main source of ionising photons, ray extended emission on scales fro m a few hundred 

although the contribution by collisionally ionised plasma to about 2 kpc has been discovered (|Wilson et al. 2000[ 

also seems to be required to accura tely match the spectra [Young et al. 200 1[ IBianchi et al. 2006|1 . X-rays bear a re- 



dViegas-Aldrovandi fc Contini 1989| ). A correlation be- mar kable morphological similarity to images in the 

tween radio power and line widths ([Wilson fc Wi llis 1980[ [[QUI] band. Simple photoionisation models applied 

IWhittle 198H |Whittle 19921 |Ulvestad fc Ho 2001[) and the to the optical to X-ray surface-brightness ratio indi- 

good m orphological correlation between [OH I] and radio cate s that solutions exist in terms of a single phase 

images ( jCapetti et al. 1996[ |Axon et al. 1998[ see however AGN-photoionised medium, where the ionisation pa- 

rameter remains constant across the whole ENLR 

Send offprint requests to: M.Guainazzi (Bi anchi et al. 2006j) . That implies a radial decrease of the 
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density cx r 2 . This is generally consistent with results 
derived from optical spectroscopy (|Kraemcr ct al. 2000, 



Table 1. Log of the XMM-Newton observations of 
NGC 1365 presented in this paper. 



Bradley et al. 2004, Krac mer et al. 2 008). although recent 



high-resolution spatially-resolved spectroscopy in a small 
sampl e of local Seyfert gala xies suggests a shallower depen- 
dence (jBennert et al. 2006P , in agreement with older results 
on Mkn 573 flCapetti et al. 1996[ ). 

High-resolution spectroscopy in the soft X-ray band 
(0.2-2 keV) confirms the overall picture. In this energy 
band, spectra of obscured AGN are typically dominated by 
recombination lines from He- and H-like transitions from 
C to Si, and by a Fe-L transitions "forest" , with negligible 
continuum contribution ([Guainazzi fc Bianchi 200"7| GB07 
hereafter). In a few bright sources, the quality of the spec- 
trum warrants detailed spectral diagnostics. In all these 
cases, photoionisation seems to be the dominant ionisation 
mechanism (Kinkhabwa la et al. 20 02, Sako ct al. 2000, 
ISambruna et al. 200T1 lArmentrout et al. 2007|) . In 
NGC 1068, by far the brightest obscured AGN of the soft X- 
ray sky, the contribution from an optically thin collisionally 
ionised plasma can be c onstrained to be ^ 1 0% of inte- 
grated soft X-ray flux (|Brinkman ct al . 2002[) . Resonant 
scattering plays also an important role in the overall ionisa- 
tion balance fKinkhabwala ~et" al. 2002| ); this constrains the 
gas column density to N H ~ 10 17 ~ 18 /Z o cm" 2 (GB07), 
where Zo is the oxygen abundance. 

In their spectroscopic study of a sample of 69 Seyfert 2 
galaxies observed with the XMM-Newton reflection grating 
system (RGS; Iden Herder et al. 200l[) . GB07 suggest that 
the conclusions derived from detailed spectral analysis of 
the brightest objects can be extended to the underlying 
population of nearby Seyfert 2 galaxies up to a flux level ~ 
10~ 13 erg cm~ 2 s -1 . However, GB07 noticed that peculiar 
objects also exist, whose X-ray spectrum exhibits different 
phcnomcnological properties from the rest of the sample. 
NGC 1365 is the X-ray brightest of these "oddities" . 

NGC 1365 is an intriguingly complex galaxy. 
Optically, it has been classified as a Seyfert either of 
type 1.5 (fVeronet al. 19801 |Hjelm fe Lindblad 1996), 1. 8 
(jMaiolino fc Rieke 1995|) . or 2~0 ([Turner et al. 1993|) . 
Following GB07 and Bianchi et al. (2009b), we 
consider in this paper NGC 1365 a "type 2" 
AGN on the basis of the column 
ing its hard X-ray emission (> 1.5 
IRisaliti et al. 20071 IRisaliti et al. 2009a|) . 
active nucl eus is highly variable at 
lengths 



cover- 

2. 



IRisaliti et a 



(Vcron et al. 1980, 



density 
x 10 23 cm 

Its obscured 
different wave- 
|Hjelm fc Lindblad 1996] 



2007). It hosts a circumnuclear star- forming 
ring with a diameter of ~1.3 kpc (~14"; 1" ~ 90 pc 
at the distance of NG C 1365: D = 18.6 ± 1.9 Mpc, 
Silb ermann et al. 1999|) . This ring is resolved into 
many compact super star clusters ([Kristen et al. 199"7[) . 
Some of them are extremely massive and young 
(jGalliano et al. 20051 IGalliano et al. 2008p. The latte r 
were detected also in the CO band (Sa kamoto et al. 2007[) . 
leading to an estimate of molecular material mass of the 
order of 10 9 Af© within 2 kpc from the AGN. The nuclear 
starburst ring must be a copious so urce of X-rays, pred omi- 
nantly produced by tw o processes {M ass-Hesse ct al. 2008, 
IGrimm et al. 20031 |Persic fe Raphaeli 2002] ). Softer 
thermal-plasma emission comes from gas heated to 
~ 10 6 K by supernova-driven shocks propagating through 
denser interstellar medium. Harder, usually featureless 



Observation # 


Start date 


Exposure time 
(ks, RGS1) 


0205590301 


17-Jan-2004 


59.0 


0205590401 


24-Jul-2004 


67.8 


0505140201 


30-Jun-2007 


126.9 


0505140401 


02-Jul-2007 


127.1 


0505140501 


04-Jul-2007 


125.1 



and asymmetrical 
IKristen et al. 1997 



X-ray emission comes from accreting compact objects 
(high-mass X-ray binaries). 

At the same time, NGC 1365 hosts a 5" bi-conical 
O Hi] outflow dHjelm fc Lindblad 1996[ 
IVeilleux et a l. 2003) and a nuclear ra- 
dio jet, the latter embedded in t he cone and extendin g 
5" SE along the galaxy major axis (Sandqvist ct al. 1995). 
The fact that the [O in] A5007 versus Ha ratio in the cone 
is larger than one suggests photoionisation by the AGN. [N 
n]/H a ratios >1 are measured out to 5-6 kpc from the nu- 
cleus (|Veilleux et al. 2003) and are consistent with either 
AGN photoionisation or shock ionisation. 

For all the above reasons, NGC 1365 represents a unique 
laboratory to study the connection between nuclear activity 
and star formation. In this paper we primarily present and 
discuss the second deepest high-resolution X-ray spectrum 
of an obscured AGN ever (after NGC 1068). It was collected 
during an occultation experiment of the NGC 1365 AGN 
in 2007 (Risal itTet al. 2009bp . Archival data from an earlier 
(2004) observation of the same field are also presented in 
this paper. 

The paper is organised as follows: we describe the obser- 
vations and data reduction in Sect. 2; in Sect. 3 we present 
a phenomenological analysis of the total RGS spectrum, 
aimed at characterising the emission line content. In Sect. 4 
we compare the RGS spectrum with physical model of opti- 
cally thin, collisionally ionised plasma and photoionisation, 
and determine the contribution of each mechanism to the 
emission line spectrum. In Sect. 5 we compare our find- 
ings with spatially resolved CCD-resolution spectroscopy 
obtained with the Chandra ACIS, and discussed in Wang et 
al. (2009). We put our results in the context of the CIELO 
(Catalogue of Ionised Emission Lines in Obscured AGN) 
sample of Seyfert galaxies (GB07) in Sect. 6, and summarise 
our main results in Sect. 7. Unless otherwise specified, sta- 
tistical uncertainties are at the 90% confidence level for one 
interesting parameter. 



2. Observations and data reduction 

The field around NGC 1365 has been observed several times 
by XMM-Newton. In order to ensure an accurate recon- 
struction of the wavelength scale in the RGS spectra, only 
observations aiming at the galaxy nucleus (Tab. [I) are con- 
sidered in this paper. 

The RGS spectra were reduced following the guidelines 
in GB07. We used the data reduction pipeline rgsproc in 
SASv7.1 (Gabriel et al. 2003)) . coupled with the most up- 
dated calibration files available in August 2008. We choose 
a fixed celestial reference point for the attitude solution 
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Fig. 1. RGS aperture (box) compared to the EPIC-pn 0.3- 
2 keV field of view (Obs. #0505 140201 in this example) 

coincident with the NED optical nucleus of NGC 1365 
(a 20 oo = 03 h 33 m 36.4s, 6 20 oo = -36 d 08 m 25s). The RGS 
aperture encompasses the whole NGC 1365 X-ray nucleus, 
including the 5 kp c extended emission imaged by Chandra 
( |Wang et al. 20 09), as well as several point-like sources in 
the galaxy (Fig. |T|) . Intervals of flaring particle background 
were filtered out by applying a fixed threshold of 2 counts 
per second on the background events, CCD#9, At = 10 s 
light curve. 

Source spectra were extracted in regions of the dis- 
persion versus cross-dispersion and Pulse Invariant versus 
cross-dispersion planes, corresponding to 95% of the point 
spread function (PSF) in the cross-dispersion direction. 
Background spectra have been generated using a sub-set 
of blank field observations, whose background count rate 
matches the level measured during each individual RGS 
observation. This method is appropriate for sources with 
extended emission on scales comparable to the PSF such 
as NGC 1365. 

The spectrum is dominated by emission lines, as typi- 
cally observed in obscured AGN (GB07) , although an con- 
tinuum emission is also clearly visible, especially in the 
wavelength range between 10 and 20A. Such a continuum 
component cannot be associated to transmission of the 
nuclear emission through a Compton-thin state, because 
the column density in these states is always higher than 
1.5 x 10 23 cm" 2 (jRisaliti et al. 20071 iRisaliti et al. 2009a[) . 
This component represents a first remarkable difference 
with respect to the archetypal Seyfert 2 galaxy NGC 1068 
( Kinkhabwala ~eT al. 2002[ ) . We compared the intensities of 
the eight brightest emission lines in the spectra accumu- 
lated during each individual observation. In all cases, the 
intensities were consistent within the statistical uncertain- 
ties. We therefore accumulated all the RGS data into a sin- 
gle spectrum ( "RGS merged spectrum" hereafter) . Its total 
exposure time is 501 and 439 ks in the RGS1 and RGS2, 
respectively. In Fig. [2]we show the spectrum obtained merg- 
ing all the RGS observations of NGC 1365. 

3. The half Ms RGS spectrum: the line content 

As first analysis step, we have characterised the emission 
line content of the spectrum. For this purpose, we fit seg- 
ments of 100 spectral channels with a simple power-law 



continuum (spectral index, L fixed to 1; the quality of the 
fit is not significantly improved if this parameter is left free 
in the fit) corrected by photoelectric absorption due to neu- 
tral gas in our Galaxy (1.3 x 10 20 cm -2 according to the 
LAB survey: Kalbcrl a et al. 2005|) and as many emission 
lines as required according to the following criterion: that 
the addition of a new emission line yields a decrease in 
the value of the C-statistics (Cash) of 9.21, formally cor- 
respondi ng to 99% confidence level for two interesting pa- 
rameter ( |Lampton et al. 1976[ ). We fit He- like triplet com- 
ponents together, imposing a difference among the compo- 
nent's centroid energy as dictated by the atomic physics. 
Lines are assumed in this paper to be unresolved, unless 
otherwise specified. 

The results are summarised in Tab. [2J where we show 
the best-fit line energy and intensity, the intensity nor- 
malised to the intensity of the forbidden component of 
the He-a Ovu ( "normalised intensity" hereafter) , the nor- 
malised intensity in NGC 1068 according to our re-analysis 
of the RGS data of this target, and the most likely iden- 
tification with the corresponding laboratory energy. The 
errors in this Table are purely statistical at the lcr level. 
Only lines detected according to the aforementioned statis- 
tical criterion are shown, except for the He-a triplets, for 
which all the components are shown, even if not all of them 
are formally detected. 

The spectrum of NGC 1365 is remarkably different from 
that of the archetypal Seyfert 2 Galaxy NGC1068. Once the 
line intensities are normalised to the forbidden component 
of the Ovu He-a triplet, we find that: 

1. Ly-a transitions of H-like species are stronger in 
NGC 1365 (except Cvi Ly-a) 

2. Fe-L transitions are stronger in NGC 1365 

3. the intensities of the He-like magnesium and silicon 
triplet components are higher in NGC 1365, although 
the poor statistics of these features in the NGC 1068 
spectrum prevents us from saying whether there is a 
difference also in the relative weight of the triplet com- 
ponents 

We'll discuss here briefly these outcomes, because they 
provide us with a guideline on the spectral analysis devel- 
oped in Sect. 4 and Sect. 5. A more detailed discussion of 
the astrophysical implications of the soft X-ray spectrum 
of NGC 1365 is, however, deferred to Sect. 6. 

The forbidden transitions of He- like triplets are strongly 
suppressed in collisionally ionised plasmas. This provides 
an explanation for Item#l above in terms of a larger con- 
tribution of a collisionally ionised plasma in NGC 1365. 
Indeed, a high ratio between the intensity of the Ly-a Ovin 
against the forbidden component of the Ovu He-a is an 
empirical criterion to distinguish starburst from AGN in 
high-resolution soft X-ray spectra (GB07). A strong con- 
tribution by Fe-L transition may also be expected in ther- 
mal plasmas (Mcwe et al. 1985). This may provide an ex- 
planation for Item#2. Nonetheless, comparatively strong 
iron lines could also be the signature of metal overabun- 
dance. This could also explain the high intensities of mag- 
nesium and silicon triplet components. Wang et al. (2009) 
found regions of high Si- and Fe-abundances N and NW 
the NGC 1365 nucleus in their analysis of the deep high- 
resolution Chandra/ ACIS image (together with regions of 
low Si- abundances primarily to the South the nucleus). We 
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Fig. 2. RGS net background-subtracted fluxed spectrum obtained after merging together all spectra accumulated during 
the XMM-Newton observations. A triangular smoothing with a 5-channel kernel has been applied for visual purposes 
only (in analogy to Fig.l in GB07). Atomic transitions observed in the RGS spectrum of NGC 1365 are labelled for 
reference. 



will later discuss the issue of heavy elements metallicity in 
the NGC 1365 interstellar medium. On the other hand, the 
intensities of the He-like nitrogen, oxygen and neon He- 
a triplet components are comparable in NGC 1068 and 
NGC 1365. 

A feature at centroid energy E ~ 738.1 eV (A = 
16.800 A) is close to the Ovn Radiative Recombination 
Continuum (RRC). This feature is narrow. If fit with a 
Gaussian profile the upper limit on its intrinsic width is 
0.8 eV. Narrow RRC features are unambiguous signatures 
of photoionised plasmas (Lied ahTet al. 1995)) . Nonetheless, 
contamination by the nearby Fexvn 2 S 2 P feature at 
16.780A is possible. We will further discuss the identifica- 
tion of this feature in the framework of the global thermal 
fit of the RGS spectrum (Sect. 4.1). 

4. Spectral fitting of the RGS spectrum 

The next step for the understanding of the physical pro- 
cesses responsible for the soft X-ray emission in NGC 1365 
is to use physical emission models to fit globally the whole 
soft spectrum or, at least, to simultaneously fit the emis- 
sion lines in Tab. [2] We have followed two complementary 
approaches: global fit of the RGS spectrum with a com- 
bination of thermal components; and simultaneous fits of 
emission line ratios through CLOUDY-based photoionisa- 
tion models. 



4.1. Thermal fit of the overall RGS spectrum 

Firstly we have tried to fit the combined RGS spectra 
with a combination of optically thin, collisionally ionised 
plasma ("thermal scenario" hereafter). We used trough this 
paper the ape c model (jSmith et al. 2001[) in Xspec vl2 
(|Arnaud 1996|) to model such a component. We used as a 
criterion to increase the model complexity that any fur- 
ther component yields an improvement of the C-statistic 
AC < —9.21. Following the results from the spatially- 
resolved spectroscopy of the Chandra diffuse emission 
dWang et al~2 009), we have assumed that each continuum 
spectral component is modified by its own photoelectric 
absorption column density. In all the models we also in- 
cluded one power-law and one blackbody continuum. These 
components are required to fit emission excesses in the 
softest (E < 0.4 keV) and hardest (E > 1.7 keV) RGS 
band. The blackbody component may account for the inte- 
grated emission of point sources (X-ray binaries and Ultra 
Luminous X-ray sources, ULX). Alternative parametrisa- 
tion for this "soft excess" are possible, such as emission 
from a disk blackbody or bremsstrahlung. They do not sub- 
stantially affect the main results present in this paper. The 
power-law component is instead required to fit a feature- 
less hard excess above ~ 1.7 keV. Only an upper limit on 
the spectral index is obtained, indicating a very hard spec- 
trum (T < —2.5). Other parametrisation of this compo- 
nent are also possible: a bremsstrahlung yields, for instance, 
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Energy (eV) 
(eV) 



Intensity (/) 

(io- 5 



i/no vim 

NGC 1365 



I/I[0 vii(f)j 
NGC 1068 



Transition (Ei a t) 
(eV) 



0.7 ±0.4 

< 0.6 

1.8 ±0.3 

0.6±§:§ 
_o. 

-0.18 

1.2 ±0.2 
1.0 ±0.2 
0.6 ±0.2 

2.4 ±0.3 

2.5 ±0.2 

i.89±8:J| 
o.47±8:il 

0.75±gi| 

.0.16 
-0.17 

o.42±8:i§ 

0.42±g;il 



0.60±!j;iH 
0.30±E!;i? 
< 0.25 

o.73±8:il 
o.27±8:i? 

0.16 ±0.07 

o.52±8:?8 
o.4i±8:?8 

0.23±8;?g 
1.00 ±0.13 

i.03±8:81 

0.78 ± 0.08 
0.20 ± 0.06 



o.3i±8:8l 

0.62 ±0.07 
0.18 ±0.05 

0.18±Hs 7 

< 0.14 

0.17 ±0.05 

o.43±8:?8 
o.24±8:8? 

< 0.25 
0.14±8;8? 

o.2i±8:8I 

< 33 

< 1.0 

0-8=1=8:1 

0.46 ±0.18 



1.16±o^ 

o.35±8:?? 

0.08 ± 0.02 

o.89±8:i4 

< 0.06 

o.n±8# 
o.4i±8:S8 
o.38±8:Sg 

0-12±8;?g 

i.oo±|j : | 

0.52±o; 9 

o.i4±8:8S 

0.104 ±0.014 



o.09±g:^ 

0.12±8;° 7 

o.06±8:8| 
o.io±8:8S 

< 0.05 

o.io±8:8J 
o.n±8:SS 

< 0.05 

< 0.04 

< 0.04 

< 0.04 

< 0.33 

< 0.05 

< 0.04 

< 0.05 



Contribution by 
thermal plasma (%) b 



367.4 ± 0.2 
430.72±8:gJ 
(r) - 4.4 
(r) - 10.8 
435.26±8ig 
459.0±8;| 
499.96±8:32 
573.99±8:8I 
(r) - 5.3 
(r) - 13.0 
653.4 ± 0.6 
726.00±8i| 
738.1±8;| 

773.54±8;§1 
824.57 ± 0.07 
872.0 ± 0.7 
921.63±8:?l 



(r) 



7.2 
13.0 



1021.13±o;?8 
1353.2±?;i 



8.0 
21.0 



(r) 

w- 

1473.0±?? 

1896.9 ±0.7 
(r) - 13.8 
(r) - 27.8 

1892.7 ± 0.8 f 



1.49± 



< 0.33 

o.40±8:^ 

1.0 ±0.2 

o.57±8:ig 

< 0.6 

o.33±8:il 
o.5o±8:il 

< 47 

< 2.4 

1.9±o:t 

1.1 ±0.4 



CVI Ly-a (367.56) 
NVI He-a (r) (430.75) 
NVI He-a (i) (426.37) 
NVI He-a (f) (419.86) 
CVI Ly-/3 (435.91) 
CVI Ly-7 (459.42) 
NVII Ly-a (500.31) 
OVII He-a (r) (574.02) 
OVII He-a (i) (569.07) 
OVII He-a (f) (561.06) 
OVIII Ly-a (653.70) 
FeXVII 3 S 2 P (726.29) 
Ovu RRC C (739.11) 

Ovm Ly-P (774.71) d 
FeXVIII 3 d 2 p (824.72) 
OVIII RRC (871.52) e 
NelX He-a (r) (922.14) 
NelX He-a (i) (915.13) 
NelXHe-a (f) (905.11) 
NeX Ly-a (1021.92) 
MgXI He-a (r) (1352.24) 
MgXI He-a (i) (1343.44) 
MgXI He-a (f) (1331.90) 

MgXII Ly-a (1472.51) 
SiXIII He-a (r) (1867.47) 
SiXIII He-a (i) (1853.61) 
SiXIII He-a (f) (1839.76) 



100 ± 50 
90 ±50 

< 100 
26 ±7 
70 ±40 
90 ±50 
100 ± 30 
110 ±30 
34 ± 12 
40 ±8 

110 ±20 
86 ± 11 
60 ±20 

83 ± 16 
98 ± 14 
60±Ig 
130 ± 60 

< 100 
90 ±30 
120 ± 30 
110 ± 15 

< 100 
90 ±50 
90 ±30 

< 100 

< 100 
90 ±30 



"most likely identification. The laboratory energy is shown in brackets, 6 contribution to the total intensity due to the best-fit 
model in Tab. 3. Values higher than 100% mean that the model over-predicts the data at that energy; c fit with a single Gaussian 
(a < 8 eV). Possible blending with Fe xvn 2 S 2 P (738.98 eV): see Sect. 4.1 for the results of a fit with two Gaussians; d possible 
blending with Fexvm 3 S 2 P (773.02 eV); e fit with a single Gaussian (a < 3 eV). Possible blending with Fexvm 2 S 2 P (872.75 eV). 
See Sect. 4.1 for the results of a fir with two Gaussians; ' unidentified. 



kT > 28 keV. Although in principle this spectral compo- 
nent could be associated with the integrated emission of 
unresolved X-ray binaries in the RGS aperture, it is diffi- 
cult to infer any meaningful physical constraints given the 
small energy range on which it yields a significant contri- 
bution to the fit. We will therefore refrain from discussing 
it any further in this paper. We only observe that the same 
component with a similar very flat spectrum is required by 
the fit of the NGC 1365 ACIS spectrum, by contrast to the 
soft continuum component which is required by the RGS 
data only (cf. Sect. 5). Finally, we included a photoelectric 
absorption component covering the whole RGS model, to 
take into account obscuration by gas in the Galaxy along 
the line-of-sight to NGC 1365. 

The model outlined above, although capable of repro- 
ducing most of the emission line intensities, fails at repro- 
ducing the Ovu and Nvn triplet ratios. This is illustrated 
in Fig. [3] which shows the spectral region around the for- 
mer feature superposed with a model with three thermal 
plasma components. As expected for the temperatures (in 
the range 0.2-1 keV) which these fits yield, the model has 
a much higher recombination component than required by 
the data. Although it is impossible to associate an absolute 



confidence level of the fit quality when the C-statistics is 
used, its high value (13042.7/5304 degrees of freedom) is 
indicative of an imperfect model. 

We have therefore removed from the fit wavelength 
ranges around the Ovu (21. 0-23. OA) and the Nvi (28.0- 
30.0 A) triplets, and repeated the fit procedure. As a cross 
check we have verified that the best-fit model obtained on 
this restricted wavelength range does not overproduce the 
observed counts in the excised wavelength ranges. 

Using the restricted wavelength range, the best-fit re- 
quires only two thermal components. In principle the abun- 
dances of each component with respect to the solar values 
(| Anders &: Grevesse 1989| have been left free to vary inde- 
pendently. However, for only a sub-set of them the corre- 
sponding confidence interval are significantly constrained. 
They are listed in Tab. [3] together with the other best- 
fit parameters. The measured Galactic column density 
{Nh,Go.i = (8±5) x 10 19 cm~ 2 ) is marginally consistent with 
the contribution estimated from the radio LAB survey. The 
C-statistic value is 6237.0/4440 degrees of freedom. In Fig. [4] 
we show the RGS spectrum and best-fit thermal model as in 
Tab. [3] superposed. The contribution to the total intensity 
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Fig. 4. RGS spectrum (crosses) and best-fit model (lines) for the whole RGS sensitive bandpass when the best-fit 
thermal model is applied. The wavelength scale is in the observer's frame. 



of each emission line due to this model ("baseline model" 
hereafter) is given in the rightmost column of Tab. O 

A multiple-temperature structure could be the results 
of fitting a single temperature non-equilibrium ionisa- 
tion plasma. We have tested this hypothesis by alter- 
natively fitting the RGS spectrum with a single vnei 
(jBorkowski et al. 2001]) model in Xspec. The resulting C- 
statistics is worse (6857.7/4422 dof). Again, this in itself is 
not enough to rule out this model. Nonetheless, we observe 
that several lines in the short wavelength regime (A & 14 A) 
are not well accounted for as in the baseline model with two 
equilibrium thermal components. Moreover, it remains true 
also with the non-equilibrium model that most of the for- 
bidden component of the Ovu and Nvi triplets are not 



accounted for. The discrepancy in the case of the non- 
equilibrium collisional scenario is even higher than in our 
baseline model. The recombination (forbidden) component 
of the Ovu He-a in the model is twice as strong (four times 
weaker) than in the data. We will not consider this solution 
any longer. 



We discuss now the identification of the feature at 
E = 738.1±g:§ cV (A = 16.799±8:So9 A )- Formally, its 
wavelength is inconsistent with the strongest close tran- 
sitions: Ovu RRC and Fe xvn 2 S 2*> (at A=16.777 and 
16.780 A, respectively). 60 ± 20% (1-cr error) of its inten- 
sity are accounted by the thermal scenario best-fit model, 
although some excess at lower energies remain. We have 
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21.5 



22 

Wavelength (A) 



22.5 



Fig. 3. RGS spectrum around the Ovu He-a triplet 
(crosses) and best-fit model (lines) in a scenario where 
only optically thin collisionally ionised plasma components 
(three in total) are used. The observer's frame wavelength 
of the forbidden (f), intercombination (i) and resonant (r) 
components is marked by the dashed lines. The apparent 
absorption feature at A ~21.8Ais due to a bad detector 
pixel. 



this additional component are: (Joviirrc = 14±g eV, and 
Ioviirrc = (9±4) x 10~ 6 photons cm~ 2 s -1 . We will refer 
to these values hereafter when discussing the properties of 
the Ovu RRC in NGC 1365. 

Similarly, the identification of the emission line at E ~ 
872 eV (A = 14.218 ± 0.011 A) is not fully unambigu- 
ous. The centroid energy is consistent with either Ovin 
RRC (A =14.228 A) of Fexvm 2 S 2P (A =14.220 A). A 
fit with two Gaussians having centroid energy fixed to 
the laboratory values, however, yields only an upper limit 
(< 4 x 10~ 6 cm" 2 s" 1 ) on the intensity of the latter. With 
this modelling the best-fit parameters of the Ovin RRC are: 
E = 870 ±2 eV; a < 8 eV; / = (1.8=bf;§) x 10" 6 cm" 2 s" 1 . 

We have also tried to constrain the bulk velocity of the 
lines, by allowing the redshift of the thermal components 
to vary. Indeed, this yields an improvement in the quality 
of the fit by AC/Ais = 121.3/1 dof for a (formal) velocity 
shift of At; = 181 ±7 km s _1 . This is comparable to the sys- 
tematic uncertainties in the RGS wavelength scale (20 mA 
at 3cr, corresponding to ~300 km s" 1 at the wavelength of 
the Ovu He-a) . We consider therefore this piece of evidence 
inconclusive as to whether the soft X-ray lines in NGC 1365 
are affected by a velocity shift of astrophysical nature. 

Assuming the baseline model as in Tab. [31 the observed 
0.3-2.0 keV flux is (2.61 ± 0.03) x 10~ 12 erg cm" 2 s _1 , cor- 
responding to a luminosity (corrected for Galactic absorp- 
tion) of 1.73 x 10 41 erg s _1 . 



tried to fit this excess with an additional Gaussian emis- 
sion line, keeping its centroid energy frozen to the source 
frame laboratory wavelength of the Ovu RRC. The fits im- 
proves by AC/Az/ = 52.3/2 dof. The width and intensity of 



Table 3. Left column: Parameters for the RGS 
Chandra/ ACIS thermal model best fits. 



and 



Parameter 



RGS 



ACIS 



Thermal components 






Nh,i (10 20 cm" 2 ) 


< 0.4 


< 1.2 


kTi (eV) 


301 ±7 


400 ± 130 


Ni a 


6.9±J;E| 


i.3±l:S 


Nh.2 (10 20 cm" 2 ) 


34±1 


< 100 


kT 2 (eV) 


690±|8 


650±g° 


N 2 a 


8.6±8:i 


i.8±U: 3 7 



Elemental abundances with respect to solar 


Zc 


o.o3i±8;8i8 


< 1.08 


Zn 


0.072±B;8H 


< 0.15 


Zo 


o.oi23±E!:S 


o.o3i±8:8?§ 


ZNe 


0.030±3;88I 


o.n±8:8S 


Z Mg 


0.083±[!;8} 2 


0.14 ±0.03 


Zsi 


o.o83±8;E!ii 


o.i5±8:88 


Zpe 


0.023 ± 0.002 


o.o43±8:8ol 


Blackbody 






kT (eV) 


364 ± 2 




N b 


3.i±8:§ 





a in units of , ,1° , — ttt f n e nndV , where Da is the angular size 
of the source and n e and uh are the electron and H densities, 
respectively 

6 in units of W~ 4 j^-, where L39 is the source luminosity in units 

10 

and D10 is the distance to the source in units of 



of 10 39 erg s" 
10 kpc 



4.2. CLOUDY-simulations of the line emission intensity ratios 

In order to estimate the contribution of photoionisation to 
the soft X-ray lines, we have generated a grid of models of 
a photoionised nebula with Cloudy (jFerland et al. 1 998). 
We have assumed a standard AGN continuum, with the fol- 
lowing parameter: "blue bump" temp erature 10 s K; g ox = 
— 1.5; ctjjv = —0.5; and a x = —1.35 (Risali ti et al. 2 000). 
where the a are spectral energy index in the band or 
between the bands specified by the subscripts. a ox was 
determined from the strictly simultaneous optical-to-X- 
ray Spectral Energy Distribution measured by the XMM- 
Newton EPIC and OM instruments during Compton-thin 
states (see Risaliti et al. 2009a, 2009b for an analysis of the 
EPIC spectra). We have generated grid of reflected spec- 
tra from a photoionised nebula, spanning a large range in 
ionisation parameter: log(t/)=[-2.0,2.0j3; electronic density 
log(n e ) = [3,14] and total column density log(A r ff ) = [17,24] 
in steps of 0.1 dex. We have extracted intensity ratios 
from these simulations for all the lines detected in CIELO 
(GB07). For each model of the grid we have calculated a 
figure of merit of the agreement between the observed and 
the simulated line ratios through a generalised x 2 function: 



E 



(h 



In 



where I Q is the observed line intensity (with its statistical 
error cr/J and Ic is the intensity predicted by Cloudy, 
both normalised to the value of the forbidden compo- 
nent of Ovu He-a. The sum is carried out on all the 
line species detected in the RGS spectrum of NGC 1365. 



the ionisation parameter U is defined as 



f(g) 



where $(i/) 



n(H)c ' 

is the surface flux of ionising photons, and n(H) is the total 
hydrogen density 
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We consider acceptable all models, which yield a x 2 value 
Y 2 < V 2 . +9.21. 

No solution exists, which can simultaneously fit the 
whole line spectrum observed in NGC 1365 in terms of one 
photoionised component, even if allowance is made for the 
elemental abundances to vary in the range [0,10 Z@\. Even 
more constraining, no solution is able to simultaneously fit 
the intensity of the Ovin Ly-a and of the recombination 
component of the Ovu He-a triplet. Fits with a combina- 
tion of two pure photoionisation models yields still xl > 5. 
On the other hand, the x 2 value corresponding to the fit 
with the RGS baseline model (i.e. two thermal components 
without any photoionisation components) is 47.0 for 19 dof, 
significantly better than for two photoionisation compo- 
nents, although still unacceptably high. A fit with one ther- 
mal and one photoionisation component is only marginally 
better: x 2 =41.5 for 19 d.o.f. (for \ogU ~1.4). 

These results, coupled with the spati ally resolved spec- 
troscopy of the Chandra NGC 1365 field flWang et al."2009| ) 
leads us to investigate in this section a physical scenario, 
where the baseline thermal model developed in Sect. 4.1 
applies, and any photoionised component is required only 
to explain the deviations between the experimental data 
and this baseline model. 

We have therefore applied the grid of Cloudy mod- 
els to a line spectrum obtained by the RGS observation 
after subtracting the contribution of the thermal scenario 
best-fit according to the rightmost column in Tab. O Only 
four of the lines detected in the RGS spectra have still in- 
tensity measurements statistically different from zero after 
the "thermal" contribution is subtracted: the forbidden and 
inter-combination component of Ovu He-a, the forbidden 
components of Nvi He-a, and (marginally) the Ovin Ly- 
j3 transition, as well as the Ovu and Ovin RRC, which, 
however, are not simulated by Cloudy. These four lines 
obviously yield only very loose constraints on the physical 
properties of the photoionised plasma responsible for their 
emission. The minimum y 2 value is 20.3 for 20 dof, corre- 
sponding to log(?7) = 1.625. In principle values of \og(U) 
comprised between 0.6 and 1.9 yield all y 2 value in the ac- 
ceptance range. One may use the contribution to the x 2 of 
the Nvi forbidden line (Fig. to further constraint, at a 
lower statistical significance, the range of acceptable ionisa- 
tion parameters. At la: log(f/) ~ 1.6±o;| (Fig- El - I n Fig. [5] 
we show iso-x 2 contour plots in the n e versus Njj plane 
for log(C7) =1.2, 1.625, and 1.925. The electron density is 
only loosely constrained to be ^ 10 10 cm~ 3 and the column 
density > 10 22 cm~ 2 at the 99% confidence level for two 
interesting parameters (these intervals are only marginally 
more tightly constrained if a less demanding confidence in- 
terval is chosen). 

We stress that these results need not to be over- 
interpreted, given the simplicity of the model assumed and 
the strong assumptions underlying it. Our results are to be 
interpreted more as a simple "proof of existence" of pho- 
toionised material in the nuclear region of NGC 1365 rather 
than as an accurate measure of the photoionised spectrum. 

5. Comparison with the Chandra results 

Deep Chandra observations of NGC 1365 
( |Wang et al. 2 009) have shown that the soft X-ray 
morphology is rather complex, with four components con- 
tributing to the overall energy output: a) the unresolved 
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Fig. 5. x 2 contribution for the Nvi (f) and Ovu (i) lines as 
a function of \og(U). 



nucleus; b) diffuse emission on scale as large as 5 kpc; 
c) off-nucleus point sources; d) the jet. Comparing the 
RGS spectroscopic measurement with CCD resolution 
spatially-resolved spectroscopy is therefore far from trivial, 
given the large RGS aperture (cf . Fig. [T]) . Only a fraction of 
this field is imaged by Chandra. Nonetheless, a comparison 
with the Chandra spatially-resolved spectroscopic results 
may be interesting to evaluate the systematic uncertainties 
induced by spatially integrated spectroscopy on such a 
complex source. We remind that the Point Spread Function 
of the XMM-Newton optics is too large for the diffuse soft 
X-ray emission associated to the circumnuclear starburst 
to be discernible. 

We have extracted a Chandra-ACIS spectrum from the 
whole NGC 1365 X-ray surface, including nucleus and point 
sources (details on the data analysis reduction can be found 
in Wang et al. 2009), and applied the best-fit RGS base- 
line model in the 0.3-2 keV energy band. The model pro- 
vides a good fit to the data (x 2 = 62.7/70 dof), and 
residuals are visible mainly at an energy consistent with 
the Ovu He-a triplet (Fig. [7]), as expected. The cross- 
calibration between RGS and the Chandra ACIS is dis- 
cussed by Plucinsky et al. (2008). Flux measurements at 
the Ovu energies by ACIS are typically lower (by ~10%) 
than the RGS. The difference becomes smaller (±5%) at 
the Neix and Nex energies. The best-fit model parameters 
are shown in the rightmost column of Tab. [3J Interestingly 
enough, the Chandra fit does not require the blackbody 
component. Taking into account the larger statistical un- 
certainties associated with the determination of the param- 
eters from the Chandra spectrum, the main differences to 
the RGS best-fit are the normalisations of the collisionally 
ionised plasma components and the metallicities. The best- 
fit metallicities measured by Chandra are systematically 
higher by a factor 2 to 4 than the RGS one (although con- 
sistent within lcr). On the other hand, the normalisations 
of the continua measured by the ACIS are around a factor 
of 5 lower than measured by RGS. We interpret these find- 
ings as due to the contribution of unresolved point sources 
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Fig. 7. ACIS total NGC 1365 spectrum (upper panel) and 
residuals against the best-fit RGS model (lower panel; best- 
fit parameters in the rightmost column of Tab. [3]) 
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Fig. 6. Iso-% 2 contour plots in the n e versus Njj plane 
for log(<7)=1.2 (upper panel) , log(C/)=1.625 (middle panel), 
and log(t/) = 1.925 (lower panel) The lines indicate the 
confidence level at 68%, 90% and 99% for two interesting 
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Fig. 8. Spectral ratio between the total NGC 1365 ACIS 
spectrum and the the ACIS spectrum obtained excising 1" 
around the nucleus. 



in the RGS aperture, which dilute the contrast between 
the emission lines and their underlying continuum. Indeed, 
point sources detected by Chandra (among them two 
ULXs, IStrateva fc Komossa 20081 ISoria et al. 2009)) con- 
tribute around two thirds of the overall integrated soft X- 
ray flux. In principle, one could sum up the spectra of these 
sources, and subtract their contribution to the RGS spec- 
trum. However, given the fact that the XMM-Newton and 
Chandra observations are not simultaneous, and that these 
sources are highly variable (jSoria et al. 2009[) . the outcome 
of this exercise would be merely academic. 

An important additional information instrumental to 
the interpretation of the RGS spectra comes from ACIS 
spatially-resolved spectroscopy. In Fig. \E\ we compare the 
best-fit total ACIS spectrum with a spectrum extracted 
from the same area once a 1" region around the nucleus 
is excised from the spectral accumulation. The main devia- 
tions occur once again at the wavelength of the Ovn He-a 
triplet and at short wavelengths. The former evidence sug- 
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Table 4. Best-fit spectral parameters when a simple pho- 
toelectrically absorbed collisionally ionised plasma model 
is applied to the NGC 1365 ACIS "total" (left column), 
and "total with AGN excised" ( "No AGN" ; right column) 
spectra 





"Total" 


"No AGN" 


N H (10 2U cm" 2 ) 
kT (eV) 
Flux a 
X 2 / dof 


9A±»i 
620 ± 20 
6.2 ± 1.3 
77.5/86 


8.8 ± 1.0 
609±Jg 

5.7 ±1.1 
72.6/73 



a in the 0.3-2 keV band in units of 10" 13 erg cm" 2 s" 1 

gests the the bulk of the photoionised plasma traced by the 
He-a emission lines is located within the nucleus unresolved 
by Chandra. The latter piece of evidence suggests that the 
overall temperature of the collisionally ionised diffuse emis- 
sion increases approaching to the centre. However, formal 
fits with a photoelectrically absorbed apec component to 
the two spectra confirms only marginally this hint (Tab. [5]). 
Only <,10% of the total soft X-ray emission is associated in 
this galaxy to the unresolved nucleus, by contrast to what 
typically observed o n AGN-dominated Seyfert 2 galaxies 
(|Bianchi et al. 2006|) . 

6. Discussion 

NGC 1365 is the brightest Seyfert 2 known, whose soft X- 
ray spectrum is energetically dominated by emission from a 
collisionally ionised plasma. This component is most likely 
associated to intense episodes of star formation occurring in 
its nuclear environment. The analysis of the RGS spectrum 
presented in this paper has provided the first proof of the 
statement above, while at the same time discovering the sig- 
natures of photoionised plasma. NGC 1365 is therefore also 
in X-rays a privileged laboratory to study the connection 
between nuclear activity and nuclear star formation. 

In this section, we put the above results in the overall 
context of the soft X-ray emission of Seyfert 2 galaxies. 

6.1. NGC 1365: the paradigm of soft X-ray weak obscured 
AGN 

In none of the four good quality soft X-ray high-resolution 
spectra of nearby AG1>Q the contribution of collisionally 
ionised plasma is significant. Whenever good quality 
spectra can be measured, they appear dominated by pho- 
toionisation, most likely due to the AGN, whose primary 
emission illuminates the NLR gas through absorption-free 
lines-of-sight (jKinkhabwala et al. 20021 ISako et al. 20001 
ISambruna et al. 20011 lArmentrout et al. 2007)) . In 
NGC 1068, by far the brightest obscured AGN of the soft 
X-ray sky, the contribution of a collisionally ionised plasma 
is constrained to be ^ 10% of the integrated 0.3-2 keV 
flux (jBrinkman et al. 2002|) . 

Arc the soft X-ray spectra of bright Seyfert 2 galaxies 
representative of the whole parent population of nearby ob- 
scured AGN? GB07 showed that the distribution of line ra- 
tios in their sample of 69 Seyfert 2s observed by the XMM- 

2 NGC 1068, Mkn 3, the Circinus Galaxy and NGC 4151 
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Fig. 9. Ratio between the Ovin Ly-a and the forbidden 
component of the Ovu He-a as a function of the integrated 
luminosity in oxygen lines for the Seyfert 2 galaxies of the 
CIELO sample (GB07). The dashed blue area indicate the 
locus occupied by a control sample of starburst galaxies 
extracted from Wu et al. (2002). Filled (Empty) circles in- 
dicate measurements (upper/lower limits). NGC 1365 is the 
filled circle with the lowest Lq value. 



Newton RGS is fairly uniform with only a few exceptions, 
NGC 1365 being the most significant of them. However, 
the statistical quality of line measurements in each individ- 
ual RGS spectrum is rather low, except for a few lines in 
a few sources. Nonetheless, GB07 discovered an empirical 
criterion, which allows us to separate on a statistical ba- 
sis the soft X-ray spectra of the CIELO Seyfert 2 galaxies 
from those of a control sample of starburst galaxies: the 
latter occupy a specific region in a plane where the inte- 
grated luminosity in oxygen lines and the ratio between the 
intensity of the Ovin Ly-a and the forbidden component 
of the Ovu He-a are plotted (Fig. 6 in their paper). We 
present a slightly clearer version of that Figure in Fig. [51 
based on a re-analysis of 92 Seyfert 2 galaxies with the 
same method and procedure described in GB07. Most of 
CIELO objects are characterised by high oxygen luminosi- 
ties (L 5 x 10 40 erg s -1 ) and low Ovin versus Ovu ratio 
(< 0.5). The starburst control sample of Wu et al. (2002) 
nicely occupies a complementary region in this plot. 

There are a few Seyfert 2 outliers, i.e. objects charac- 
terised by low oxygen luminosities and high Ovin versus 
Ovu ratios. They are: NGC 1365, NGC 4303, NGC 5506, 
and NGC 7582. Apart from NGC 1365, NGC 4303 and 
NGC 7582 are well known to host intense episode of star 
formation close to the nucleus (Bianchi et al. 20071 and ref- 
erences therein) , which in NGC 4303 dominate the UV and 
soft X-ray energy output CJimcncz-Bailo n et al. 2 003) . It is 
therefore tempting to identify these outliers with sources, 
whose soft X-ray spectrum is substantially "contaminated" 
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Table 5. Emission line intensities in stacked RGS spectra 
of oxygen luminosity selected, low counts CIELO sources. 



Transition 


Low oxygen-L a 


High oxygen- 


Fe xvm 3 s 2 P 


0.47 ±0.26 


< 0.14 


Fe xvn 3 d 2 3 


0.67 ±0.34 


< 0.20 


Fe xvn 3 s 2 P 


1.00 ±0.50 


< 0.23 


Fe xxi 3 d 2 P 


0.76 ± 0.43 


< 0.95 


Ne IX He-a (r) 


0.57 ±0.32 


0.17 ±0.10 


Ne IX He-a (i) 


< 0.20 


< 0.09 


Ne IX He-a (f) 


0.24 ±0.22 


0.21 ±0.10 


N VI He-a (r) 


< 0.85 


< 0.30 


N vi He-a (i) 


0.70 ± 0.47 


0.20 ±0.17 


N vi He-a (f) 


0.93 ±0.56 


0.35 ±0.18 


N vil Ly-a 


0.72 ± 0.44 


< 0.17 


Vin Ly-a 


1.25 ±0.58 


0.48 ±0.14 


vin Ly-/3 


0.44 ±0.26 


0.48 ±0.14 


vil He-a (r) 


0.84 ±0.50 


0.52 ± 0.20 


vil He-a (i) 


< 0.54 


< 0.24 


vil Hc-a (f) 


1.00 ±0.41 


1.00 ±0.20 



a normalised to the Ovu He-a (f) 



ative weakness of the AGN photoionisation contribution 
in these objects, as well as the overall weakness of their 
soft X-ray emission lines, are due to shielding of the ion- 
ising high-energy AGN emission by optically thick matter 
in the BLR. In this respect, however, NGC 5506 does not 
seem to match this scenario. Large (by a factor of about 
2) historical variation of the X-ray flux are accompanied 
by a very modest variability of the absorbing column den- 
sity (Guainazzi et al., in preparation). The column den- 
sity of NGC 5506 (Bianchi et al. 2003j) is also significantly 
smaller than measured even in the Compton-thinnest states 
of NGC 1365 and NGC 7582 (N H ~2.5xl0 22 cm" 2 ). It is 
still possible that such a small column density is associated 
to matter in the host galaxy rather than in the nuclear envi- 
ronment (jPiconcelh et al. 20071 IMatt et al. 2003|) , and that 
the ENLR gas is looking at the AGN trough a different opti- 
cal path where obscuration of the active nucleus is substan- 
tially high. As already pointed out by several authors, we 
are only now starting to glimpse the true complexity of the 
gas and dust distribution in the innermost parsec around 
supermassive black holes (Elvis 2000, Risalit i et al. 20 05b, 
lElitzur fe Shloshman 2006j) ~ 



by collisionally ionised optically thin emission associated to 
star formation. 

One may try and make another step forward by adding 
together the information from all those spectra, whose indi- 
vidual data quality is not enough to warrant significant line 
detection. We do that by stacking together rest frame spec- 
tra of CIELO sources with total background counts in the 
range between 50 and 500. In Tab.[5]we show the normalised 
intensity of the emission lines detected in the stacked spec- 
tra for objects with total integrated oxygen lines luminosity 
smaller and higher then Lo ~ 10 41 erg s _1 . Errors in this 
Table are at 1-a. Low-oxygen luminosity sources exhibit: 
a) stronger Fe-L lines; b) stronger recombination and inter- 
combination transitions in He-a triplets of Nexi and Nvi; c) 
stronger Ly-a transitions of Ovin and NeVH All the above 
trends are consistent with a stronger role played by colli- 
sionally ionised plasmas in low-oxygen luminosity sources. 

6.2. BLR shielding of the ionising continuum? 

CIELO Seyfert 2 which are under-luminous in oxy- 
gen lines are not necessarily intrinsically low-luminous. 
If we consider those galaxies, which have a measure- 
ment of the Ovin Ly-a versus Ovn(f) intensity ratio in 
Fig. [5J and exclude NGC 4303, where the presence of 
an AGN is still under debate ( Jimencz-Bailon et al. 2003): 
the absorption-corrected 2-10 keV luminosity is 2.5 x 
10 42 erg s- 1 in NGC 1365 (Risaliti et al. 2009aP; 1.5- 
3.2xl0 42 erg s" 1 in NGC 7582 (jBianchi et a l. 2009a), and 
4.5-9.1 xlO 42 erg s" 1 in NGC 5506 (jBianchi et al. 20031) . 
NGC 1365 and NGC 7582 share another remarkable obser- 
vational property: both exhibit extreme (by one order of 
magnitude) variability of the X-ray absorbing column den- 
sity on timescales as short as 10 hours ([Risaliti et a l. 2009b) 
and ^ 1 day (jBianchi et al. 2009ap . respectively. Such rapid 
transitions from Compton-thin (transmission-dominated) 
to Compton-thick (reflection-dominated) states can be ex- 
plained only if the X-ray absorbers are located at distances 
from the black hole comparable to that of the broad line 
region (BLR). It is intriguing to speculate that the rel- 



6.3. Properties of the photoionised plasma in NGC 1365 

It is impossible to significantly constraint the properties of 
the photoionised plasma responsible for the bulk of the for- 
bidden component of He-like triplets in the RGS spectrum 
of NGC 1365. With the assumed SED: log(tf) = 1.6±g-| 
(Ict), n e <, 10 10 cm" 3 , and N H ^ 10 22 c m - 2 . From the 
best-fit CLOUDY solution (Sect. 4.2), and using the def- 
inition of ionisation parameter, we estimate that the in- 
ner side of the photoionised nebula is constrained to be 
located <; 0.75 pc from the source of the ionising photons. 
This is still consistent with a location of this gas within 
the BLR. However, it is not consistent with the soft X- 
ray photoionised gas being the same responsible for the 
highly ionised, outflowing iron absorption lines discovered 
by XMM-Newton (jRisaliti et al. 2005a|) . and probably as- 
sociated to a disk outflow. The width of the oxygen RRC 
features constrains the gas temperature to be ~10 5 K, of 
the same order of magnitude as the temperature measured 
in NGC 1068 (jKinkhabwala et al. 2002) or other CIELO 
AGN (GB07). The only He- like ion, for whose He-a triplet 
meaningful constraints on the diagnostic parameters G and 
R ( jPorquet fe Dubau 2000| can be derived is Ovu. The lat- 
ter is primarily sensitive to density. R = f /i = 3.6± 1.2 
implies log(n e ) <, 10.5, in agreement with the results of the 
comparison between Cloudy simulations and the RGS line 
spectrum. On the other hand, G = (/ + i)/r > 4.2 trans- 
lates into a loose constraints on the photoionised plasma 
temperature, ^ 9 x 10 5 K, again in agreement with the 
RRC-based diagnostic. 

6.4. Properties of the collisionally ionised plasma in 
NGC 1365 

A formally good fit to the RGS spectrum of NGC 1365 re- 
quires two optically thin, collisionally ionised thermal com- 
ponents with temperatures ~300 eV and ~640 eV. These 
values should not be interpreted too literally, though. The 
best-fit temperature is mainly driven by the fit of the emis- 
sion line strength and position. The tightest constraints are 
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Fig. 10. 6.2/im PAH feature flux (in arbitrary units) versus 
the 2-10 keV absorption corrected luminosity in a sample 
of obscured (grey) and unobscured (black) AGN. NGC 1365 
is the dot surrounded by an empty circle. 



provided by the strongest lines, which are mostly emitted 
at ~0.3 keV (e.g: oxygen and neon lines) and ~0.7 keV 
(Mg lines). Gas at various other temperatures may be also 
present, however it could not get the same weight in the 
fitting process because it does not emit as strongly in those 
lines. It is likely that this results should be interpreted as a 
distribution of temperature over the nominal best-fit range. 

Our analysis qualitatively confirms the results by Wang 
et al. (2009) that the gas in the circumnuclear environment 
of NGC 1365 has an overall very low metallicity. However, 
RGS lacks spatial resolution matching the complex mor- 
phology of metal abundances in the NGC 1365 nuclear am- 
bient gas. Besides the overall metal under-abundance, re- 
gions of silicon and iron over-abundance were discovered by 
Chandra ( |Wang et al. 2 009). Moreover, the high Equivalent 
Width of the relativistically broadened iron K a line mea- 
sured by the XMM-Newton EPIC (jRisaliti et al. 2009a|) is 
also consistent with a factor of about three iron overabun- 
dance. 

Gas in the spiral arm of the galaxy should be have an 
even lower metallicity. At the same time, NGC 1365 has a 
high number of ULX in its disk (jSoria et al. 2009]) . These 
pieces of evidence are consistent with a scenario, whereby 
ULXs are black hole remnants of very metal-poor (Z < 
O.IZq) stars, which little mass in their winds (m ~ Z ' 7 ). 

Does NGC 1365 host a more luminous starburst when 
compared with obscured AGN of comparable intrinsic AGN 
power? In Fig. [TU] we plot the flux of the PAH feature 
at 6.2/x as measured by Spitzer as a function of the un- 
obscured 2-10 keV luminosity as measured by CAIXA 
(jBianchi et al. 2009b[) in a sample of obscured and unob- 
scured local Universe AGN. PAH features are a strong 
indicator of star formation in active galaxies (see, e.g., 
ILutz et al. 20081 and references therein). NGC 1365 is the 



second brightest PAH AGN in the sample after NGC 7469; 
the latter is well known to exhibit an extremely bright star- 
burst at distances between tens to about 10 3 of parsecs from 
the AGN (jDi'az-Santos et al. 2007) . 

7. Conclusions and summary 

In this paper we primarily discuss the 0.5 Ms RGS spectrum 
of the Seyfert 2 galaxy NGC 1365. The main results of our 
analysis can be summarised as follows: 

— no variability either in the overall spectral shape or in 
the line intensity is detected when comparing observa- 
tions distant up to five years apart 

— the spectrum is dominated by emission lines from op- 
tically thin, collisionally ionised plasma, plus a strong 
continuum primarily in the 10-20 A range. Formally, 
a good fit is obtained with two components, with 
temperatures ~300 eV and ~640 eV, the latter seen 
through a photoelectric absorption column density ~ 
3 x 10 21 cm" 2 . These results are in broad agreement with 
the range of temperatures and column densities mea- 
sured by spatially resolved moderate resolution spec- 
troscopy of this galaxy with Chandra ( Wang et al. 2009) 

— the RGS spectrum confirms the suggestion by Chandra 
that on the average sub-solar metallicities are required 
to fit the data, although the unknown level of contami- 
nation by individual point-like sources makes the exact 
measurements of metallicities from the RGS spectrum 
uncertain 

— a residual contribution by a photoionised plasma is re- 
quired to fit the forbidden component of He-like tran- 
sitions of oxygen and nitrogen, as well as the inter- 
combination component of Ovu. The physical param- 
eters of the photoionised plasma can be only loosely 
constrained: log(J7) ~1.6, n e <J 10 10 cm -3 , and Njj *S 
10 22 cm~ 2 , and correspond to a plasma at a distance 
<; 0.75 pc. The detection of oxygen RRC features sug- 
gests a photoionised temperature ~10 5 K, in agree- 
ment with values typically observed in AGN-dominated 
Seyfert 2 galaxies. 

Despite its comparably low luminosity in soft X-ray lines 
with respect to AGN-photoionisation dominated Seyfert 2 
galaxies, NGC 1365 is not intrinsically X-ray weak. We 
speculate that the soft X-ray spectrum NGC 1365 could ap- 
pear as starburst-dominated due to shielding of the primary 
AGN ionising continuum by optically thick matter associ- 
ated to the BLR, a scenario which may apply to other AGN 
(such as NGC 7582), which exhibit extreme and rapid vari- 
ability of the absorber column density. Good-quality soft 
X-ray spectra of "changing-look" AGN would be crucial to 
test this hypothesis. Alternatively, NGC 1365 could host 
a remarkably luminous nuclear starburst when compared 
to its AGN accretion power, as suggested by the luminous 
PAH features detected in its Spitzer spectra. 
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